ON THE PERIODIC SOLUTIONS OF A RIGID 
DUMBBELL SATELLITE IN A CIRCULAR ORBIT 
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■^^ , Abstract. The aim of the present paper is to provide sufficient con- 

ditions for the existence of periodic solutions of the perturbed attitude 
dynamics of a rigid dumbbell satellite in a circular orbit. 

1. Introduction and statement of the main results 

Q \ In this paper we consider the attitude dynamics, perturbed by small 

(H torques, of a rigid dumbbell satellite in a circular orbit under the gravita- 

tional torque of a central Newtonian force field. Our objective is to provide 
sufficient conditions for the existence of periodic motions about the satellite's 
center of mass that are asymptotic to translational motion in an absolute 
coordinate system. This type of motion, denoted, cylindrical equilibrium, 
is well known in the astrophysics literature on satellite's dynamics (see for 
instance [HEHl]). 

^sO \ These motions have an important application to satellite orientation prob- 

lems because a satellite can reach some specified nominal regime along pe- 
riodic trajectories only through the influence of gravitational torques and 
other small perturbed torques induced by some control mechanism. 

m ' Following the methods developed in [8] , the equations of motion governing 

the attitude dynamics of a rigid dumbbell satellite are 
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with 6 and (/> the eulerian angles of nutation and precession. The perturbed 
torques F* , are smooth functions periodic in the variable t with 
, / dd dd>\ , / dQ dd>\ 

and e a small real parameter. In this work we are interesed in the periodic 
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functions emerging from the equilibrium solution 9 = and </> = of ([TJ 
when e — > 0. 

By means of the change of coordinates x = 8 and y = 4> and linearizing 
the equations of ([!]) in this equilibrium we obtain 

(2) 
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and fj smooth functions in the variables I t, —r, —r I • On the other hand 
Jl y'dt'dtj 

the functions /j are periodic in t and in resonance p:q with some of the 
periodic solutions of the unperturbed dumbbell satellite given by 

d 2 x 

— + 3x = 0, 
(3) dt 

d y 

The objective of this work is to provide, using as a main tool the averaging 
theory, a system of nonlinear equations whose simple zeros provide periodic 
solutions of the perturbed dumbbell satellite in circular orbit with equations 
of motion given by ([2]) . Some other works using similar techniques (see Llibre 
et al. [T] for more details), are [2[ [3] or [4] where the method is used for action 
angle variables. In order to present our results we need some preliminary 
definitions and notations. 

The unperturbed system ([3]) has a unique singular point, the origin with 
eigenvalues 

±y/3i, ±2i. 

Consequently this system in the phase space I x, ~r,y, -j~ ) has two 

planes filled with periodic solutions with the exception of the origin. These 
periodic solutions have periods 

nn 2vr m 

T 1 = - or T 2 =vr, 

according to whether they belong to the plane associated to the eigenvectors 
with eigenvalues ±v3i or ±2i, respectively. We shall study which of these 
periodic solutions persist for the perturbed system ([2]) when the parameter 
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e is sufficiently small and the perturbed functions F{ for i = 1, 2 have period 
either pTi/q, or pT2/q, where p and q are positive integers relatively prime. 
We define the functions: 

•Fi(X 0) Yo) = T- [ P1 sin (V3i) A x (t) /i (t,A 2 (t),0)dt, 
&W Jo v / 

F 2 (X ,Y ) = y— coslVst) A^t) f 1 (t,A 2 (t),0)dt, 

Gi(Z ,W ) = — / sin(2t) A 3 (i) f 4 (t,0,A 4 (t))dt, 
PK Jo 

2 f" pT ' 2 
g 2 (Z ,W ) = — / cos(2i) A 3 (t) / 4 (i,0,A 4 (t))dt, 

with 

Ai(t) = X cos (V3t) + -p sin (>/3t) , 

v3 

A 2 (t) = Y cos (\/3i) - y^Xo sin (>/3t) , 
A 3 (i) = Z cos (2t) + -^ sin (2t) , 

A 4 (t) = Wo cos (2t) - 2Z sin (2t) . 
A zero (Aq, Y *) of the nonlinear system 
(5) J- 1 (A o ,Y ) = 0, F 2 (X ,Y ) = 0, 

such that 

' d{F u F 2 ) 
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(A- ,yo)=(A*,Y *)/ 

is called a simple zero of system ([5]). Similarly, a zero (Zq,Wq) of the 
nonlinear system 



(6) Gi(Z o ,W o ) = 0, g 2 (z ,w ) = o, 

such that 

' d(g 1 ,g 2 ) 



det 



5(Z ,Wo 



^0, 
(^o,Wb)=(^J,WJ)/ 

is called a simple zero of system ([6]). 

Our main results on the periodic solutions of the perturbed dumbbell 
satellite ([2]) are the following. 

Theorem 1. Assume that the function Ff and F 2 of the perturbed dumb- 
bell satellite with equations of motion |7P are periodic in t of period pT\ / q 
with p and q positive integers relatively prime. Then, for e 7^ suffi- 
ciently small and for every simple zero (Xq,Y q *) ^ (0,0) of the nonlin- 
ear system ([5]), the perturbed dumbbell satellite |2P has a periodic solution 
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0{t, e), — (t, s),4>(t, e), — (t, e) J with 

lim (e(Q,e),^(0,e),<l>(Q,6),^(0,e)\ = (X*,Y*,0,0). 

Theorem [I] is proved in section [2j Its proof is based in the averaging 
theory for computing periodic solutions, see the Appendix of the paper. 

We provide an application of Theorem [JJ in the following corollary, which 
will be proved in section [3j 

Corollary 2. We consider the system $T\) with 

d0 , dd>\ „ , y-v . AQ „ fdO\ 2 



^(M,^^fJ=c OS ^-sin(V3t)sin^--sin^-j 
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Then, for e/ sufficiently small the system (QP /ias one periodic solution 
with 

lim ^(0, £ ),|(0, e ),^(0,,),^(0,,)) = (f ,0,0,o) . 

Similarly we obtain the following result. 

Theorem 3. Assume that the functions F£ and F£ of the perturbed dumb- 
bell satellite with equations of motion |7]) are periodic in t of period pT<ijq 
with p and q positive integers, relatively prime. Then for e ^ suffi- 
ciently small and for every simple zero (Zq,Wq) ^ (0,0) of the nonlin- 
ear system (|6|) ; the perturbed dumbbell satellite (Q]j has a periodic solution 

0(t, e), — (t, e),4>(t, e), — (t, e) J with 

lim ^(O,e),^(O,e),0(O,e),^(O,e)J = (0,0, Z* ,W£). 

The proof of Theorem [3] is analogous Theorem [TJ 

On the other hand, we provide an application of Theorem [3] in the next 
corollary, which will be proved in section [3j 

Corollary 4. We consider the system (QP with 
„ „( n d9 , d(f>\ , . n d4 , , n . , ( d(j)\ d<p 

/ n d9 . d<j)\ . ± . ._ . M . ' ' N :! 
2 \ ~dt ~dt) =sm ^" sm ( 2t ) sin ^" sin Vi r// 
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Then, for e/0 sufficiently small, the system f7]J has three periodic solutions 
with 

Jim («(O, e ),f(O, e ),0(O, e ),|(O, e) ) = (o,0,l,^), 
Jim (e(0,e), f t (0,e), 0(0,8), ^(0,e)^j = Uo,l^E,o\ , 

Jim ^(0,e),^(0, e ),^(0,e),^(0,e)) = U,0,i±^,0J . 

Remark 1. The momenta apply to a rigid dumbbell satellite in a circular 
orbit under the gravitational torque of a central Newtonian force field are, 
in general, functions of the Euler's angle and some maps which depend of 
the independent variable t (time). The applications presented in Corollaries 
\M and \jj\ represent some usual momenta which model the solar radiation 
obtained by the solar panels located at the satellite, see for more information 

2. Proof of Theorems Q] and \3\ 

(dx dv \ 
x i—r:yi~r\ we write the 
at at J 

differential system of the perturbed dumbbell satellite (J2]) as a first-order 

differential system defined in R 4 , as follows: 



dX 
~~d~t 


= Y, 


dY 




lit 




dZ 
~dt ' 


= W, 



3X + eF l (t, X, Y, Z, W) + e 2 Ri (t, X, Y, Z, W, e) 

(7) 



^- = -AZ + eF 2 (t, X, Y, Z, W) + e 2 R 2 (t, X, Y, Z, W, e) . 
at 

System ([7]) with e = is equivalent to the unperturbed dumbbell satellite 
system ([3]). On the other hand, the periodic orbits of the unperturbed 
system are described in the following lemma. 

Lemma 5. The periodic solutions of the unperturbed system with e = are 

X(t) = X cos (y/3t) + -p sin (y/3t) , 

v3 

( 8 ) Y(t) = Y cos (V3t) - V3X sin (\/3i) , 

Z{t) = 0, 
W(t) = 0, 
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of period T\ , and 



(9) 



of period T 2 ■ 



X{t) 
Y{t) 

Z{t) 



0, 

o, 



Wo 



Z cos (2i) + ^p sin (2t) , 



W(t) = Wo cos (2t) - 2Z sin (2t) 



Proof. Since the unperturbed system is a linear differential system, the proof 
is routine. □ 

Proof of Theorem d Assume that the function F\ and F 2 of the dumbbell 
satellite system with equations of motion $Z§ are periodic in t of period pTi/q 
with p and q positive integers, relatively prime. Then system (|2j) is periodic 
in t with period pT\. 

We shall apply Theorem [6] of the appendix to the differential system ([7j). 
We note that system (|7|) can be written as system (lllj) . taking 
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I ° \ 

F 1 (t,X,Y,Z,W) 



\ F 2 (t,X,Y,Z,W) J 



and 



G 2 (t,x,e) 



/ \ 

e 2 i?i (t, X, Y, Z, W, e) 



V e 2 R 2 (t,X,Y,Z,W,e) J 

We shall study which periodic solutions (jSJ) of the unperturbed system 
corresponding to the system ([7]) with e = can be continued to periodic 
solutions of the unperturbed system for e/0 sufficiently small. 

We shall describe the different elements which appear in the statement of 
Theorem [6J Thus we have that Q = M 4 , k = 2 and n = 4. Let n > be 
arbitrarily small and let r 2 > be arbitrarily large. We take the open and 
bounded subset V of the plane Z = W = as 



V = {(X , Y , 0, 0)eK 4 :n< ^X 2 + Y 2 < r 2 } 
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As usual C1(V) denotes the closure of V. If a = (Xq,Yo), then we can 
identify V with the set 

{a£l 2 :ri< ||a|| < r 2 }, 

here || • || denotes the Euclidean norm of R 2 . The function /3 : C1(V) — > M? 
is /3(a) = (0, 0). Therefore, in our case the set 



Z = {z a = (a, (3(a)) , a G C\(V)} = {(X , Y , 0, 0)eK 4 : n < ^X 2 + Y 2 < r 2 }. 

Clearly for each z a £ Z we can consider the periodic solution x(£,z a ) = 
(X(t), Y(t), 0, 0) given by © of period P T X . 

Computing the fundamental matrix M Za (t) of the linear differential sys- 
tem with e = associated to the T-periodic solution z a = (Xq, Yq, 0, 0) such 
that M Za (0) is the identity of M 4 , we conclude that M(t) = M Za (t) is equal 
to 



/ 



cos (v3i) 



sin 



(V3t) 



\ 




-V^sin (\/St) cos(V3i) 

cos(2£) 








sin(2t) 



\ 



-2 sin (2 1) cos(2£) / 

Note that the matrix M Za (t) does not depend of the particular periodic 
solution x(t,z a ). Since the matrix 



/ 




M _1 (0) - M~ 1 (pT 1 ) 








l-cos(^f^) 



sm 






/ 4y / 3p7T 



\ 



\0 



sin(^ 



i_ cos r4v^r 



/ 



satisfies the assumptions of statement (ii) of Theorem [6] because the deter- 
minant 



1 - cos (^f^ 



sm 



\ 3 



sin(^ 



1_ cos (4v|pzE 



4si n 2 ( 2 ^)/0, 



we can apply this theorem to the unperturbed system. 

Now f : M 4 -> R 2 is £(X, F, Z, W) = (X, F). We calculate the function 



g(X ,Y Q ) = G(a) = i 



\pTi 



i Jo 



M-Ht)Gx(t,x(t,z a ))dt 
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and we obtain 

/ Ji(x 0j y ) 



( 1 /•^sin( v / 3t) 



PTljQ 



V^ 



Fi(t J Ai(t) J A 2 (t),0 ) 0)(ft 



i -i-/ 1 cos(V3t) J Fi(t,A 1 (i),A 2 (t),0,0)dt 



Using 

iq(t, Ai(t) ; A 2 (t), 0, 0) = Ai(t) h (t, A 2 (t), 0) 

we obtain the functions given by (j3|). Then, by Theorem [6] we have that for 
every simple zero (Xq,Yq) £ V of the system of nonlinear functions 

(10) Ji(X o ,y o ) = 0, F 2 (X ,Y ) = 0, 

we have a periodic solution (X, Y, Z, W)(£,ff) of the unperturbed system 
such that 

(X, Y, Z, W)(t, e) -+ (J^, Y *, 0, 0) as e ->■ 0. 

□ 



3. Proof of the two corollaries 

Proof of Corollary \^ Under the assumptions of Corollary [2] the nonlinear 
equations (|SJ) becomes 

Y (3X%+Y 2 ) 



Fi(X ,Y 



M*o,Y 



48 



(3( v / 3 - 3X )X$ - (v/3 + 3X )y 2 ) 



24 



This system has the following real solution 

(X o * ! F o *)=(^,0 



Moreover 



/ 



det 



d{F u F 2 ) 



d{X ,Y ) 



1 

384 



(*o.*b)=(#.°) 

check that this solution is simple. So, by Theorem Q] we only have one 
periodic solution of ([T|) . This completes the proof of the corollary. □ 

Proof of Corollary [7j Under the assumptions of Corollary [5] the nonlinear 
equations (|5J) becomes 



Gi(Z ,W 



G2(Z ,W ) 



(go - l)Wb 

8 
(4Z (2 + Zq - 2Zg) - W 2 (l + 2Z )) 
32 
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This system has the following four real solutions 
(Z* ,WZ)=U±^Y (Z* ,W£) = 
and 

The solutions which differ in a sign are different initial conditions of the 
same periodic solution of the system ([2]) and 



det 



det 



det 



d(Gi 


,9i) 


d(z 


W ) 


d(Gi 


,Qi) 


d(z 


W ) 


d(Gi 


,Qi) 


d(z 


W ) 



1 

(Z ,Wo)=(l,±^) I 32 ' 



17-7^17 

\ 17 + 7VU 
{ZoiWo)= (l±m,o) I 512 



Therefore, by Theorem [3] we only have three periodic solutions of the per- 
turbed dumbbell satellite. This completes the proof of the corollary. □ 

Appendix: Basic results on averaging theory 

In this appendix we present the basic result from the averaging theory 
that we shall need for proving the main results of this paper. 

We consider the problem of the bifurcation of T-periodic solutions from 
differential systems of the form 

(11) x(t) = G (t,x) + eGi(t,x)+e 2 G 2 (t,x,e), 

with e = to e 7^ sufficiently small. Here the functions Go , G\ : R x SI — >• 
R n and G2 : R X Vt x (— e , £0) — > R n are C 2 functions, T-periodic in the first 
variable, and f2 is an open subset of M n . The main assumption is that the 
unperturbed system 

(12) x(t) = G (t,x), 

has a submanifold of periodic solutions. 

Let x(i, z, e) be the solution of the system (fT2|) such that x(0, z, e) = z. We 
write the linearization of the unperturbed system along a periodic solution 

x(t, z, 0) as 

(13) y(i) = AcG o (£,x(i,z,0))y. 
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In what follows we denote by M x (t) some fundamental matrix of the linear 
differential system (flgjh and by £ : R fe x M. n ~ k — > R k the projection of R n 
onto its first k coordinates; i.e. £(xi, . . . ,x n ) = (%l, ■ ■ ■ ,%k)- 

We assume that there exists a ^-dimensional submanifold Z of f2 filled 
with T-periodic solutions of (|12p . Then an answer to the problem of bifur- 
cation of T-periodic solutions from the periodic solutions contained in Z for 
system (jlip is given in the following result. 

Theorem 6. Let V be an open and bounded subset of Mr, and let j3 : 
C1(V) — > W l ~ k be a C 2 function. We assume that 

(i) Z = {z a = (a, /3(a)) , a G C1(V)} C £1 and that for each z a £ Z the 

solution x(i, z Q ) of (PT21) is T-periodic; 
(ii) for each z a G Z there is a fundamental matrix M Za (t) of (|13p such 
that the matrix M - ' 1 ^) — Af~ 1 (T) has in the upper right corner the 
kx(n—k) zero matrix, and in the lower right corner a (n—k)x(n—k) 
matrix A a with det(A a ) ^ 0. 

We consider the function Q : C1(V) — > R 

(14) g( a ) = z(J-J M- a \t)Gx(t,x(t,z a ))dty 

If there exists a £ V with G(a) = and det ((dQ /da) (a)) 7^ 0, then there is 
a T-periodic solution <p(t, e) of system (|lip such that ip(0, e) — > z a as e — > 0. 

Theorem [6] goes back to Malkin [5] and Roseau [6] , for a shorter proof see 
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